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Abstract 

This paper presents a class of Two-Step General Linear Methods for the numerical 
solution of Retarded Functional Differential Equations. Explicit methods up to order five 
are constructed. To avoid order reduction for mildly stiff problems the uniform stage order 
of the methods is chosen to be close to uniform order. 



1 Two-step General Linear Methods for 
Ordinary Differential Equations 

For the numerical approximation of the solution y{t) of a system of Ordinary Differential 
Equations 

y'{t) = f{t,y), te[to,T], 

where / : R x M'^ — ^ R'^, yo e R"^, 

we consider the class of General Linear Methods [5] 

= ^ a,, + j2 "'^ ^ = 1, . . . , ^, 

i=i j=i 

= E hpf^ + E ^^^-^^r'^' ^ = 1, • • • > (2) 

i=i i=i 
Ff] = /(t„_i + Q/i,F/"]), z = l,...,s, 

where y\'' ^\ ■ ■ ■ ^ yr' — input vectors, available at step number n, 

. . . , fJ"' — stage values, f["\ . . . , fJ"' — derivative values, a^-, Uij, bij, Vij. — coefficients 
of the method. 
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Let us restrict ourselves to two-step GLMs and choose r = s + 2, 



~ y{tn-i), yt ^' ~ y{tn-2), y^"/' ~ hy'{tn-2 + Cih) , i = l,...,s. Then 

yl""^^ = yn-l, = yn-2, yfc^' = hf{tn-2 + Q/i, F-'""^') = /iF/""^^, f = 1, . . . , S. 

and (El) takes the form 

s s 
= ^ X] ^j""^ + ^iil2/n-l + Ui2yn-2 + Ui^2+j Fj''~^\ i = l,...,S, 

i=i j=i 

s s 
yn = h^ bij F]"' + Vuyn-l + t^l22/n-2 + /i ^ ^^1, 2+^^]""^' , (3) 

i=i j=i 
Ft^ = f{tn^, + cAYt^), z = l,...,s. 

In the construction of GLMs it is assumed that |/|"~^' = Uiy{tn-i) + Vi hy'{t„-i) + 0{h'^) and 
'preconsistensy conditions' holds 

:::::: ^ 

For ([3]) we have Ui = 1, Vi = 0, M2 = 1, f 2 = —1, "^2+4 = 0, V2+i = 1, i = 1, . . . , s. 
It follows from jl) that 

tii2 = l-^^il, Z = 1,...,S, 
Vi2 = l- Vu- 

Let us denote 



' := 1<I ' ttij := Ui^2+j, Ui:=Uii, ^ Ui2 = I - Ui, j = 1, . . . , s, i = 1, . . . , s, 
bj:=bij, bj:=vi^2+j, v := Vu, ^Vi2 = l~v, j = l,...,s, 

then the method ([3]) satisfying 'preconsistensy conditions' ([4]) takes the form 

s s 

y^= {l-v)yn.2 + vyn-i + hJ^bjK^r'^ + hY,hKf\ 
Kfl = /(t„.i + Q/.,F/"]), ^ = l,...,s, (5) 

s s 

y.["l = (1 - Ui)yn-2 + + /i ^ ajji^]""^^ + ^ aijKf\ i = 1, . . . , s. 



2 Two-step General Linear Methods for 

Retarded Functional Differential Equations 



We begin with notations introduced in [1]. 

Let r G [0, +oo), and C be the space of continuous functions [— r, 0] — > W^, equipped with the 

maximum (uniform) norm = max |(^(^)|, <^ ^ C, where | ■ | is an arbitrary norm on MJ^ . 

ee[-r,o] 

Let u be continuous function [a — r, b) — > M'^, where a < b. Then Vt G [a, b) shift function is 
well defined by Ut{9) = u{t + 6), 9 e [-r, 0], and Ut G C. 

Let us consider a system of Retarded Functional Differential Equations 

y'it) = fit.yt), te[to,T], 

y,^{e) = cf>{e), 9e[-r,0]. 

where (to, 4>) eVL, f : Q — > M"', Q cR x C, Q is open set. 
It is assumed that there exists a unique solution of Q. 

We introduce the class of two-step GLMs for RFDEs on the base of approach proposed 
in [1]. We can reformulate the method ^ for RFDEs (JG]) as follows 



r^M(a/i) = (1 -t;(a))r7["-^l(0) + v{a)r]^''-'\h) + /i ^ 6,(«)k]"-'^ + /i J] 6,(a)irf', 

j=i j=i 
«G [0,1], 

irf]=/(t„_i + Q/.,F[t:J, z = l,...,s, 

s s 

yM \ah) = (1 - Mi(a))r7["-il(0) + Ui{a)r]^"-'^\h) + /i ^ ai,-(a)ir]"~'^ + /i ^ ai,(a)K]"l, (7) 

a G [0,q], i = l,...,s, 

where 

77'"^^' : [— r, h] — > W^, Kf'~'^^ G are available as approximations 

computed in the step — 1, 
Y^^^ * — stage functions, 

k\^^ — stage values, 

Wi(-), aij(-), aij(-), f(-), &j(-)5 &j(-) ^ coefficients of the method. 



[n-1] [n] 



« a ^//////////////^ ^ 

tn-1-h-r t„_,-h t„_, t„_, + h 

Thus r/W(a/j) ~ + ah), ae [0, 1], r]^''^{9) ^ + ^ G [-r,0], hence 

^i"' ~ on [-r - /i, 0], where t„ = t„_i + /i. 

Remark 2.1. VFe chose two-step methods among multi-step methods (k > 2) for the following 
reasons. 

• For multistep methods, the local error E{h, t„-i, ?/t„_i) has the required order only if exact 
solution y{t) is sufficiently smooth on [tn-k, tn]- This is rather severe assumption for many 
problems ((61). For the case of two-step methods (k = 2) this condition imposes the weakest 
restriction on stepsize. 

• Furthermore, in the case of two-step methods, starting procedure and stepsize strategy 
seem to be simplest ones. 

Let us denote 

77 := r^W, ,= r) := r^I-^l, := Kt'\ a := 

then the method dT]) can be reformulated in Stefano Maset's notations as follows 



r]{ah) = {1 - v{a))r]{0) + v{a)ri{h) +h'^hj{a)Kj h^bj{a)Kj, a e [0,1], 

j=i j=i 

K, = f{a + cAYl,), 

s s 

Y\ah) = (1 - Ui{a))r]{0) + u^{a)r]{h) + h^aij{a)Kj + h^aij{a)Kj, a E [0,Ci], (8) 

Vi0) = vd0), Oe[-r,0], 

Y\e) = r],{e). ee[-r,o], 

^ K. 



, ^//////////////^ ^ 

a-h-r a- h a a + h 



3 Two-step GLMs for RFDEs in Stefano Maset's notations 



Let us consider a system of Retarded Functional Differential Equations 

x'{t) = f{t,xt), te[to,T], 



(9) 



where (to, <p) E il, f : — > M"*, C M x C, is open set. 

It is assumed that conditions of existence and uniqueness theorem for the (l9|) hold. 

When s-stage Two-Step General Linear Method for RFDEs (TSGLM) with coefficients 
(aij(-), bj{-), Ci,aij{-), bj{-),Ui{-), t'(-))ij=i,...,s is applied with stepsize /i to (I9|) for the computation 
of the solution x{t), it yields, as an approximation on [— r, h] of the shift function y := x{a + ■), 
the function 



r]{ah) = (1 - v{a))rj{0) +v{a)r]{h) + h^bj{a)Kj + h^bj{a)Kj, a G [0, 1], 
v{0)=VhiO), Oe[-r,0], 



(10) 



where 



function ij/^ ^ y on [— r — h,0] and stage values Ki, are available as approximations 
computed in the previous step, 



K, = f{a + cAY, 



1 = 1, 



fill 



and : [— r, Cj/i] — > is a stage function given by 

s s 

Y\ah) = {1 - Ui{a))f}{0) + Ui{a)f}{h) + h'^aij{a)Kj + aE [0,Ci], (12) 



It is assumed that coefficients {aij{-),bj{-),Ci,aij{-),bj{-),Ui{-),v{-))ij=i^,,,^s of TSGLMs 



satisfy the following conditions: 



o-ijV)^ O'ijV)^ 



are polynomial functions [0, c,] 

bj{-), bj{-), v{-), are polynomial functions [0, 1] 

Ci G M, Ci> 0, 

ay(0) =a,,(0) = 0, Ui{0) = l, 

6,(0) =6,(0) =0, t;(0) = l, 



h 3 = 1, 

J = l, 
1 = 1, 

h 3 = 1, 



s. 

s. 
s. 
s. 

s. 



(13) 

(14) 
(15) 
(16) 



The last two conditions correspondingly gurantee continuity of the stage functions Y^,f^ G C and 
the approximate solution rjh E C provided that approximate solution computed in the previous 
step is continuous function rjf^ E C. 



Remark 3.1. If the conditions 

Ui{-) = 1, aij{-)=0, i,j = l,...,s, (17) 

hold, the two-step method ^ becomes the one-step RK method for RFDEs introduced in [1], 
where initial function (p :=r]f^. 

Definition 3.2. TSGLM with coefficients (ajj(-), bj{-), Ci,aij{-), bj{-), '^(■))i,i=i,.--,s is called 
explicit if ajj(-) = for all j : j > i, i, j = 1, . . . , s. 

Definition 3.3. The function E = rj — y : [0,h] — > M.'^ computed under the assumption that 
fj^ = y on [-r - h, 0] is called the local error of TSGLM ((HI). 

Definition 3.4. The function E"" = — y : [0,Cj/i] — > M.'^ computed under the assumption 
that r]f^ = y on [— r — h, 0] is called the local stage error of TSGLM (l8|). 

By analogy with [7] we define stage order for RFDEs. 

Definition 3.5. Let functions i?* = — y : [0, Cih] — > R'^, i = 1, . . . , s, 

^s+i — ^ _ y . j^Q /jj — ^ j^d g^j,g computed under the assumption that 

fjf^ = y on [— r — h,0], Kj = y'{—h + Cjh), Kj = y'{cjh), j = 1, . . . , s, that is 

s s 

E\ah) = (1 — Ui{a))y{—h) + Mj(a)y(0) + /i^^ aij(a)?/'(— /i + Cjh) + h aij{a)y'{cjh) — 

-y{ah), ae[0,Ci], i = l,...,s, (18) 

s s 

E'+\ah) = {1 - v{a))y{-h) + v{a)y{0) + h^bj{a)y'{-h + Cjh) + h^bj{a)y'{cjh)- 
j=l j=i 
— y{ah), a G [0, 1]. 

Denote c^+i := 1. If there are positive integers Pi and reals Di > 0, H > such that 

max \E'{ah)\< DihP'^\ he[0,H], i = l,...,s + l, (19) 

a&[0,Ci] 

then positive integer p = min{pi, . . . ,Ps+i} is called uniform stage order of TSGLM (El). 



4 Order conditions 



Assume that / is of class C' with respect to the second argument for a sufficiently large I and 
solution x{t) of ([9]) is of piecewise class for a sufficiently large m. 

We introduce the polynomial functions Ffc : [0, 1] — > M and Tik : [0, q] — > M given by 



{k-l)\ 



{i-v{am-iy 

k 



+ ^6,(a)(-(l-c,))^-i +J2b 

«e [0,1], 



(20) 



{k-l)\ 



+ 5^ «,,(«)(-(! -c,))^-i + 5^ 
a G [0,Ci], i = l,...,s 



a 
"k 



Remark 4.1. // the conditions ( 1171 ) hold the Fjfc, are the same as for the one-step RK 
method [1]. 

Let c^, . . . , c*. such that < < • ■ • < c*, and {c\, . . . , c*.} = {ci, . . . , Cs}, i.e. c* are 
distinct q in increasing order. 

Lemma 4.2. Lei p he a positive integer. If x is of piecewise class and the local stage 

errors computed in the previous step are = 0{h^), i = 1, . . . , s, then the local error E and 
the local stage errors E'^ satisfy 



E{ah) = hJ2bj{a)D, +^|/('=)(0)/i'=Ffc(a) a G [0, 1] 



k=l 
P 



E\ah) = hJ2 + Yl y^'\0)h''T,kia) + 0{W^^), « G [0, c,], z = 1, . . . , s, 

j=i k=i 

where 

Di = f{(J + Cih, yc,h + ^cU) - /(^ + Ci^' y^^h), 
*The Lemma [42] has been proved, but its proof is omitted here for brevity.* 



i = l,...,s. 



In the following we assume that the TSGLM satisfies the conditions Fi = and 
Fji = 0, i = 1, . . . , s, that is 



f211 



f22) 



(23) 



Via] 



1 +^6,(«) +J2bj 



a) = a, 



Ui{a) - 1 + ^aij(a) + ^ 



aij[a) = a, 



ae [0,1], 



a G [0,Ci], 2 = l,...,s. 



(24) 



The above condition is an equivalent form of uniform stage order one condition. 
Theorem 4.3. A TSGLM satisfying f l24l ) has uniform order two iff V2 = 0. 

*The theorem 14.31 has been proved, but its proof is omitted here for brevity.* 

Theorem 4.4. Let TSGLM satisfy ( l24l) and has uniform order two. 

// Ts = and E f'i(a)ri2(/5) = 0, ae[0,l], Pe[0,c*J, m = l,...,s* 

then the method has uniform order three. 

*The theorem 14.41 has been proved, but its proof is omitted here for brevity.* 

Theorem 4.5. Let TSGLM satisfy ( l24l) and has uniform order three. 

If r4 = 0, 

s 

5^ 6,(«)r,3(/3) = 0, aG[o,i], /?e[o,c:;], m = l,..., 

s s 

5^ 5^6,(a)a,,(/3)r,2(7) = 0, a e [0,1], f3e[0,c*J, 7e[0,cr], l,m = l,..., 
i=l^ j=l 

'^i—'^m Cj=C* 

then the method has uniform order four. 

*The theorem 14.51 has been proved, but its proof is omitted here for brevity.* 

Theorem 4.6. TSGLM has uniform stage order p iff 
r^fc = 0, Tfc = 0, i = l,...,s, k = l,...,p. 

Proof. Follows by Taylor series expansion of functions i?* given by ( fTSl ). 

The following results can be obtained as corollary of theorems (14.61) and (14.41) . (14.51) 

Corollary 4.7. Let TSGLM has uniform stage order two. 
If Fa = then the method has uniform order three. 

Corollary 4.8. Let TSGLM has uniform stage order three. 
If F4 = then the method has uniform order four. 

The results of Corollary (14.71) and (14. Sp can be easily generalized as follows. 

Theorem 4.9. Let TSGLM has uniform stage order p. 
It has uniform order p = p+ 1 iff Fp+i = 0. 

*The theorem 14.91 has been proved, but its proof is omitted here for brevity.* 



5 Construction of explicit two-stage GLMs of 
uniform stage order four and five 



Consider two-stage explicit TSGLM satisfying (l24|) . It's Butcher tableau is 



Table 5.1. Butcher tableau for 2-stage explicit TSGLMs 



Cl 


Ui 


(a) 


an 


(a) 


ai2 


(a) 








C2 


U2 


(«) 


^21 


(«) 


022 


(a) 


^21 (a) 







v{ 


a) 


bi{ 


a) 


hi 


a) 


bi{a) 





A natural choice will be to space out the abscissae q, i = uniformly in the 

1 s — 2 

interval [0, 1] so that [6] Ci = 0, C2 = . . . , Cg-i = -, Cg = 1. In the case of s = 2 we 

have Cl = 0, C2 = 1. 

Since Ci = 0, conditions Tik{a) = 0, a G [0,ci], k = 1,2,... reduce to rij!c(0) = 0, 
k = 1,2, . . . that follows from (fTSl) . It also follows that Wi(-) = 1, aii(-) = 0, ai2(-) = 0. 

For brevity we omit the argument a of the method coefficient functions. By theorem (14.91) . 
the method has uniform order four and uniform stage order three if 
Ffc = 0, /c = 1,2,3,4 and r2fc = 0, k = 1,2,3, that is 



-(1 


-v) 


1 


— V 




2 


1 


— V 




3 


1 


— V 




4 


(1- 


- U2) 


1 - 


- U2 




2 


1 - 


- U2 



a" 



hi + 62 = y, 



'2 



'J' 



hi +h2 = — , (26) 



- ^21 

+ a2i 



a, 
a" 



2 ■ 
3 



a 
"3 



where a G [0, 1]. 



The coefficients are defined by 



U2 = -(2a- l)(a + l)% 
V = {a — if {a + 1)^ , 
«2,i = a^ (a + 1) , 

&i = -^a'(« + l)(5a-7), (27) 

0-2,1 = a (a + 1)^ — a2,2, 

1 2 ~ 

6i = — a (2 a — 3) (a + 1) — 62, 
3 

^2 = ^ (a + 1)^ , 



where 022, ^2 remain free. The relation r5(l) = — 7^ implies that it is impossible to attain 

15 

discrete order five. 

The uniform order and the uniform stage order can be increased by finding a suitable value 
for C2. Assume that ci = 0, C2 7^ (in general case C2 7^ 1). By theorem (14.91) . the method has 
uniform order five and uniform stage order four if 
Ffc = 0, A; = 1,2,3,4,5 and Tsfc = 0, A; = 1,2,3,4, that is 



-(1 -v)+bi +62 




Vhi +b2 


= a, 


\ (1 
2 ^ 


£2)^2 


+ 02^2 


~ 2 ' 


-'3%^ +(1- 


C2fh2 


+ clb2 


a' 
~ 3 ' 


^7-?. -d- 


C2)% 


+ clb2 


a^ 
~ 4 ' 


5 


C2)% 


+ c% 


a' 

~ 5 ' 


-(1 - U2) +0.21 +0.22 




\- 021 


= a, 


I-U2 ^ , 

2 «21 (1 


02)0.22 




~ 2 ' 


1 — W2 ~ 

^ + a2i + (1 - 


€2)^022 




^ y 


4 «21 (1 


C2fo.22 




~ 4 ' 



where a G [0, 1] in the first five equations (l28l) and a G [0, C2] in other ones. 



The coefficients are defined by 
,2 ( 

= (a + 1) 1 - 2 a + 



2c2-l, 

(a + 1)^ ((10 a - 5) - 15 caa^ + (a + 1) (6 _ 3 « + 1)) 

5C22-1 ^ 

2. , (a + l)^(3c2-l) 

a2i = a a + 1 — , 

2C2 (2c2 - 1) 

022 



2C2(C2-1)(2C2-1)' 

{a + 1) (20c2^- (30a + 10) cs^ + {Vla^ ^"ia-Vi) + {^a^ + \\a+Z) C2 - 2a (a + 1)) 

4C2 (5C22 - 1) (C2 + 1) ' 

a^ (a + l)^(5c2^ - (4a -3) C2 — 2 a) 



021 = a (a + 1) I 1 



4C2 (5C22-1)(C2-1) 

a (3c2-2) 



(29) 



2(2C2-1)(C2-1) 

a (a + 1)^ (20c2^- (30a + 20) c^^ + (12a2 + 21a-4) c^^ + (-4a2 + 3a+4) C2 - 2a (a + 1)) 

4C2 (5C22 - 1) (C2 - 1) 

o? (a + 1)^ (5 ca^ - (4 a + 7) C2 + 2 a + 2) 



4C2(5C22-1)(C2 + 1) 

To attain the discrete stage order five, we determine C2 from r25(l) = 0. We have 



11-V41 , , 

c. = (30) 



16 (17 -2v^) ^ 

The relation r6(l) = ; -;=t- 7^ implies that it is impossible to attain discrete 

75(71-11^41) 

order six. 

So we construct explicit TSGLM of uniform order five, uniform stage order four and discrete 
stage order five. 

Remark 5.1. There is not a method of uniform stage order two in a class of explicit one- 
step RK methods for RFDEs. Indeed, for explicit one-step RK methods Ci = 0, C2 7^ and 

«2,j = 0, j = 2,. . . ,s, hence T2k = --^ 7^ 0, a G (0, C2], A; = 2, 3, . . . . 

It is known [7] that methods with low stage order suffer from the order reduction phenomenon 
when applied to stiff ODEs. Hence, explicit TSGLMs may he more appropriate for some mildly 
stiff RFDEs ( of course, if the smoothness conditions in Remark \2.1\ and in the begining of 
section\^ hold). 
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